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Large numbers of ground states of two-dimensional Ising 
spin glasses with periodic boundary conditions in both direc- 
tions are calculated for sizes up to 40^^. A combination of a 
genetic algorithm and Cluster-Exact Approximation is used. 
For each quenched realization of the bonds up to 40 indepen- 
dent ground states are obtained. 

For the infinite system a ground-state energy of e = 
— 1.4015(3) is extrapolated. The ground-state landscape is 
investigated using a finite-size scaling analysis of the dis- 
tribution of overlaps. The mean-field picture assuming a 
complex landscape describes the situation better than the 
droplet-scaling model, where for the infinite system mainly 
two ground states exist. Strong evidence is found that the 
ground states are not organized in an ultrametric fashion in 
contrast to previous results for three-dimensional spin glasses. 



Introduction In this work two-dimensional Edwards- 
Anderson (EA) ± J spin glasses Q are investigated. They 
consist of N spins ai ~ ±1, described by the Hamiltonian 



(1) 



The spins are placed on a two-dimensional (d=2) 
square lattice of linear size L with periodic boundary 
conditions in both directions. Systems with quenched 
disorder of the nearest-neighbor interactions (bonds) are 
investigated. Their possible values are Jy = ±1 with 
equal probability. A constraint is imposed, so that 

A question which has led to much controversy is, 
whether many pure states exist for finite dimensional, i.e. 
realistic spin glasses. For the infinite ranged Shcrrington- 
Kirkpatrik (SK) Ising spin glass |^ this question was 
answered positively by the replica-symmetry-breaking 
mean-field (MF) scheme by Parisi Additionally this 
solution of the SK-model exhibits an ultrametric 

structure: The distances d"^ between the states do not 
only satisfy the triangular inequality d"^ < d'^'^ + d'^^ but 
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the stronger ultrametric inequality d"^ < max{d°''^ jd"^^) 
as well. For an introduction to ultrametricity see 
Numerical work on the subject for the SK-modcl can be 
found in @J|. 

A complete different model is proposed by the Droplet 
Scahng (DS) theory [p[-p^. It suggests that only two 
pure states (related by a global flip) exist and that the 
most relevant excitations are obtained by reversing large 
domains of spins (the droplets). Some authors ||lj] dis- 
believe that the MF theory is consistent for realistic spin 
glasses but the question about the existence of many pure 
states is not answered. 

In this work the point is not addressed by investigat- 
ing spin glasses at finite temperature, but a different ap- 
proach is used: an analysis of ground states is performed. 
The existence of many pure states implies that between 
two ground-state configurations arbitrary differences are 
possible. Otherwise two ground states would only differ 
by the spin orientations in some finite domains, which is 
always possible in the ± J model because of the discrete 
structure of the interaction distribution. Such a simple 
structure was found for example for the 3d random-field 
Ising model p5[ | . An investigation whether the ordering 
of the ground states is ultrametric is of interest on its 
own: it is conjectured that an ultrametric structure pro- 
hibits that the usual cluster algorithms for simulating a 
system at finite temperature work efficiently . 

Next a short overview about numerical results of 
finite-dimensional systems for investigating the low- 
temperature behavior are given. In four dimensions the 
MF picture is well established. Even numerical evidence 
for ultrametricity at finite temperature but below the 
transition temperature Tq was found |]l7| , |l8| . 

For three-dimensional systems recently evidences for 
the validity of some basic features of the MF picture 
were found using simulations at finite temperature p9| ] 
and ground-state calculations . First attempts to find 
ultrametricity by simulation at finite temperature are 
given in [ pT| , p2[ . Evidences for an ultrametric ground- 
state landscape of 3d systems were recently found . 

The higher-dimensional systems exhibit a spin-glass 
phase at nonzero temperature |2j,|2^, but the 2d model 
orders only at T = . So the question arises whether 
it is possible to detect this difference in the structures of 
the ground-state landscape as well. Several results us- 
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ing direct calculations of ground states of 2d spin glasses 
using approximation methods are known |2^-j3^. Also 
exact ground-states up to L = 50 have been analyzed 
1^,^. But in all these studies only one ground state 
per realization was calculated, so it was not possible to 
study the structure of degeneracy. First attempts of the 
calculation of many different ground states can be found 
in nil. 

In the present study many different and independent 
true ground states of 2d ± J spin glasses are examined in 
detail. A combination of a genetic algorithm |3^,^^ and 
Cluster-Exact Approximation (CEA) Jsst is used. Re- 
cently it has been shown that this method is able to cal- 
culate true ground states of spin glasses |^ . Since true 
ground states are calculated directly, one does not en- 
counter ergodicity problems or critical slowing down like 
in algorithms which base on Monte-Carlo methods. The 
ground-state landscape is investigated and it is discussed 
whether the MF or the DS picture describes its structure 
better. Also the ground states are analyzed whether they 
are ordered in an ultrametric fashion. 

The paper is organized as follows: at first the algo- 
rithm used here for the calculation of the ground states 
is explained. Next all observables are defined. Then the 
results are presented. In the last section a conclusion is 
driven. 



ALGORITHM 

For readers not familiar with the calculation of spin- 
glass ground states now a short introduction to the sub- 
ject and a description of the algorithm used here are 
given. A detailed overview can be found in |38 



FIG. 1. The simplest frustrated system: a triple of spins, 
each pair of spins connected by antiferromagnetic bonds 
(dashed lines). It is not possible to satisfy all bonds. 

The concept of frustration |3|] is important for under- 
standing the behavior of ± J Ising spin glasses. The sim- 
plest example of a frustrated system is a triple of spins 
where all pairs are connected by antiferromagnetic bonds, 
see fig. 1^. A bond is called satisfiedii it contributes with a 
negative value to the total energy by choosing the values 
of its adjacent spins properly. For the triangle it is not 
possible to find a spin-configuration were all bonds are 
satisfied. In general a system is frustrated if closed loops 
of bonds exists, where the product of these bond- values 
is negative. For square and cubic systems the smallest 



closed loops consist of four bonds. They are called (ele- 
mentary) plaquettes. 

As we will see later the presence of frustration makes 
the calculation of exact ground states of such systems 
computationally hard. Only for the special case of the 
two-dimensional system with periodic boundary condi- 
tions in no more than one direction and without external 
field a polynomial-time algorithm is known |po[| . Now for 
the general case three basic methods are briefly reviewed 
and the largest system sizes which can be treated are 
given for three-dimensional systems, the standard spin- 
glass model. 

The simplest method works by enumerating all 2^ 
possible states and has obviously an exponential run- 
ning time. Even a system size of 4^ is too large. The 
basic idea of the so called Branch- and- Bound algorithm 
1^ is to exclude the parts of the state space, where no 
low-lying states can be found, so that the complete low- 
energy landscape of systems of size 4'^ can be calculated 

il- 

A more sophisticated method called Branch- and- Cut 
p3| , p3| works by rewriting the quadratic energy function 
as a linear function with an additional set of inequali- 
ties which must hold for the feasible solutions. Since not 
all inequalities are known a priori the method iteratively 
solves the linear problem, looks for inequalities which are 
violated, and adds them to the set until the solution is 
found. Since the number of inequalities grows exponen- 
tially with the system size the same holds for the com- 
putation time of the algorithm. With Branch-and-Cut 
anyway small systems up to 8'^ are feasible. 

The method used here is able to calculate true ground 
states up to size 14"^. For two-dimensional systems 
sizes up to 50^ can be treated. This is about the same 
size Branch-and-Cut can solve, but in contrast to that 
method the algorithm used here is able to calculate many 
independent ground states for each realization of the ran- 
domness. The method bases on a special genetic algo- 
rithm [^,0 and on Cluster- Exact Approximation [P5| . 
CEA is an optimization method designed specially for 
spin glasses. Its basic idea is to transform the spin glass 
in a way that graph-theoretical methods can be applied, 
which work only for systems exhibiting no frustrations. 
Next a description of the genetic CEA is given. 

Genetic algorithms are biologically motivated. An op- 
timal solution is found by treating many instances of the 
problem in parallel, keeping only better instances and 
replacing bad ones by new ones (survival of the fittest). 
The genetic algorithm starts with an initial population 
of Mi randomly initialized spin configurations (= indi- 
viduals) , which are linearly arranged using an array. The 
last one is also neighbor of the first one. Then Uq x Mi 
times two neighbors from the population are taken (called 
parents) and two new configurations called offsprings are 
created. For that purpose the triadic crossover is used 
which turned out to be very efhcient for spin glasses: a 
mask is used which is a third randomly chosen (usually 
distant) member of the population with a fraction of 0.1 
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of its spins reversed. In a first step the offsprings are 
created as copies of the parents. Then those spins are se- 
lected, where the orientations of the first parent and the 
mask agree jiZj. The values of these spins are swapped 
between the two offsprings. Then a mutation with a rate 
of p„i is applied to each offspring, i.e. a fraction p,„ of 
the spins is reversed. 

Next for both offsprings the energy is reduced by ap- 
plying CEA: The method constructs iteratively and ran- 
domly a non-frustrated cluster of spins. Spins adjacent 
to many unsatisfied bonds are more likely to be added to 
the cluster. During the construction of the cluster a local 
gauge-transformation of the spin variables is applied so 
that all interactions between cluster spins become ferro- 
magnetic. 




FIG. 2. Example of the Cluster-Exact Approximation 
method. A part of a spin glass is shown. The circles repre- 
sent lattice sites/spins. Straight lines represent ferromagnetic 
bonds the jagged lines antiferromagnetic interactions. The 
top part shows the initial situation. The construction starts 
with the spin at the center. The bottom part displays the fi- 
nal stage. The spins which belong to the cluster carry a plus 
or minus sign which indicates how each spin is transformed, 
so that only ferromagnetic interactions remain inside the clus- 
ter. All other spins cannot be added to the cluster because it 
is not possible to multiply them by ±1 to make all adjacent 
bonds positive. Please note that many other combinations of 
spins can be used to build a cluster without frustration. 

Fig. ^ shows an example of how the construction of 
the cluster works using a small spin-glass system. For 2d 
± J spin glasses each cluster contains typically 70 percent 
of all spins. The non-cluster spins act like local magnetic 
fields on the cluster spins, so the ground state of the clus- 
ter is not trivial. Since the cluster has only ferromagnetic 
interactions, an energetic minimum state for its spins can 
be calculated in polynomial time by using graph theo- 
retical methods p5[-^: an equivalent network is con- 
structed iQ, the maximum flow is calculated p9| , |50| p] 
and the spins of the cluster are set to their orientations 
leading to a minimum in energy. This minimization step 
is performed nmin times for each offspring. 

Afterwards each offspring is compared with one of its 
parents. The pairs are chosen in the way that the sum 
of the phenotypic differences between them is minimal. 
The phenotypic difference is defined here as the number 
of spins where the two configurations differ. Each parent 
is replaced if its energy is not lower (i.e. not better) 
than the corresponding offspring. After this whole step is 
done Uo x Mi times, the population is halved: From each 
pair of neighbors the configuration which has the higher 
energy is eliminated. If more than 4 individuals remain 
the process is continued otherwise it is stopped and the 
best individual is taken as result of the calculation. 

The representation in fig. ^ summarizes the algorithm. 

The whole algorithm is performed times and all 
configurations which exhibit the lowest energy are stored, 
resulting in uq statistically independent ground-state 
configurations. 

This algorithm was already applied to examine the 
ground-state landscape of 3d spin glasses . 



^Implementation details: We used Tarjan's wave algorithm 
together with the heuristic speed-ups of TrafF. In the con- 
struction of the level graph we allowed not only edges {v,w) 
with level(u') = level (u) -1-1, but also all edges {v,t) where t is 
the sink. For this measure, we observed an additional speed- 
up of roughly factor 2 for the systems we calculated. 
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algorithm genetic CEA({ Jj^}, Mi, Uo, Pm, "min) 
begin 

create Mi configurations randomly 

while (Mi > 4) do 

begin 

for i = 1 to no X Ali do 

begin 

select two neighbors 

create two offsprings using triadic crossover 
do mutations with rate Pm 
for both offsprings do 
begin 

for j = 1 to riniin do 

begin 

construct unfrustrated cluster of spins 
construct equivalent network 
calculate maximum flow 
construct minimum cut 
set new orientations of cluster spins 
end 

if offspring is not worse than related parent 
then 

replace parent with offspring 

end 
end 

half population; Mi ~ Mi/ 2 
end 

return one configuration with lowest energy 
end 

FIG. 3. Genetic Cluster-exact Approximation. 



OBSERVABLES 

For a fixed realization J ~ {Jij} of the exchange in- 
teractions and two replicas {erf}, {erf}, the overlap |^ is 
defined as 

i 

The ground state of a given realization is characterized 
by the probability density Pj{q). Averaging over the 
realizations J, denoted by [■]av, results in {Z — number 
of realizations) 

J 

The probability densities are symmetric as no external 
field is applied: Pj{q) ^ Pj{-q) and P{q) = P{-q). 
Hence, only averages of are relevant: 

j\rPj{q)dq (4) 
W^j\rP{q)dq (5) 



The Droplet model predicts that only two pure states 
exist, implying that P{q) converges for L ^ oo to P{q) = 
^{^{q ^ Qea) + S{q + qEA)), while in the MF picture the 
density remains nonzero for a range —qsA ^ 9 ^ Qea 
with peaks at i^max (0 < gmax < Qea, gmax qEA for 
L oo). Consequently the variance 

a'{\q\)^ j'^iW\~\q\rmdq=W-k\' (6) 

stays finite for L ^ oo in the MF pictures ((J^(|g|) = 
Soo + si/L"^) 0, while CT^dgl) L^y ^ according 
the DS approach. Here, y is the zero-temperature scaling 
exponent which is denoted as 6 in pl| , p^ . 

Another way of describing the finite-size behavior of 
is to sum up the contributions from small overlap- 
values q < qo'- 

P{q)dq (7) 

-la 

This value should converge to in the DS picture as long 
as go < qEA while it should stay non-zero for the MF 
framework. 

The overlap defined in (j^) can be used to measure the 
distance d"^ between two states: 

= 0.5(1 - q°'P) (8) 

with < c?"'' < 1. For three replicas a,/?, 7 the usual 
triangular inequality reads d""^ < d"'^ + d"^^ . Written in 
terms of q it becomes 

q^fi > q»-y + q^'^ ~1 (9) 

In an ultrametric space Q the triangular inequality 
is replaced by a stronger one d"^ < ma,x{d°''^ ,d'^^) or 
equivalently 

g"^ > min(g"T,9T'^) (10) 

An example of an ultrametric space is the set of leaves 
of a binary tree: the distance between two leaves is de- 
fined by the number of edges on a path between the 
leaves. 

Let qi < q2 < 93 be the overlaps q"'^ , q"'' , q'^^ ordered 
according their sizes. By writing the smallest overlap on 
the left side in equation (^0|), one realizes that two of the 
overlaps must be equal and the third may be larger or 
the same: 91 = 92 < 93 

In a finite-size system this relation may be violated. 
Here two ways are used of determining whether ground 
states of realistic spin glasses become more and more 
ultrametric with increasing size L: 

• The difference 

Sq = q2-qi (11) 
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is calculated for all triplets. Because the influ- 
ence of the absolute size of the overlaps should 
be excluded the third overlap is flxed: 173 = qjix- 
In practice only overlap triples are used where 
93 G b/iKj 9/1x2] holds to obtain sufficient statis- 
tics . With increasing size L the distribution P{Sq) 
should tend to a Dirac delta function for an ultra- 
metric system 

• If two overlaps are fixed (g"''' = q^'' — q/ix, in prac- 
tice q"'^,q'^'' e [qfix,qfix2]), equation (^) implies 
q = q"^ > 2qfix — 1 while ultrametricity implies 
Q > qjix which is stronger if qfix < 1 The 
distribution P2-fix{q) of the third overlap is used 
to characterize the ultrametricity of a system. Ad- 
ditionally the weighted fraction of the distribution 
outside [qfix.qEA] 

Il^J P2-f,x{q){q~qf,xfdq (12) 

+ / P2-ftx{q){q - qEAf dq 

(see jlj]) should vanish for L 00 in an ultramet- 
ric system. 



springs per configuration, nmin = number of CEA mini- 
mization steps per offspring, r = average computer time 
per ground state on a 80MHz PPC601. 




FIG. 4. Average ground-state energy of 2d ±J Ising spin 
glass for linear dimensions 5 < L < 40 and FSS fit of form 
e{L) = 600 +eiL'^^ resulting in e^o = -1.4015(3). 



RESULTS 

We used simulation parameters determined in the fol- 
lowing way: For each system size several different combi- 
nations of the parameters Mi,no,nmin,Pm were tested. 
For the final parameter sets it is not possible to obtain 
lower energies even by using parameters where the cal- 
culation consumes four times the computational effort. 
Using parameter sets chosen this way genetic CEA cal- 
culates true ground states, as shown in p5| ]. Here a mu- 
tation rate oipm = 0.05 and nn = 40 runs per realization 
were used for all system sizes. Table 1 summarizes the 
parameters and gives the typical computer time r spent 
per ground state computation on a 80 MHz PPC601. 
For each system size ground states for 1000 different re- 
alizations of the disorder were calculated. On average 
uq > 29 (not necessarily different) ground-state configu- 
rations were obtained for all system sizes L using nu = 40 
runs per realization. 



(sec) 



5 8 
10 16 
14 16 
20 32 
32 128 
40 512 



0.02 
0.4 

3 

30 
780 
5400 



Tab 1. Simulation parameters: L = system size, Mi 
= initial size of population, Uo — average number of off- 



At first the result for the ground-state energy as func- 
tion of system size is presented in fig. ^. By perform- 
ing a finite-size scaling (FSS) analysis using the function 
e(L) = Coo + eiL~'^'^ a ground-state energy of the infinite 
system of Coo = —1.4015(3) is obtained. This is consis- 
tent with an extrapolation from exact ground states of fi- 
nite systems Coo = —1.4015(8) and with other former 
results from transfer-matrix calculations too — —1.402(1) 
Monte-Carlo simulations = -1-407(8) ||8|, mul- 
ticanonical simulations eoo = —1.394(7) [pO|, genetic al- 
gorithms Coo = -1.400(5) (31I, Coo = -1.401(1) ||La 
special cluster-construction method Coo — —1.402(2) |2g|] 
and pure CEA eoo = -1.400(5) |3|]. The value presented 
here has a higher accuracy than the former results. 

The exponent of the decay is 62 — 2.2(1) which is much 
larger than the upper limit oi y = 62 < {d — l)/2 ~ 0.5 
proposed by DS. This is an evidence that the DS picture 
may not be appropriate for 2d realistic spin glasses. 

Information about the ground-state landscape can be 
extracted by evaluating the distribution of overlaps. In 
fig. H f'(lgl) is shown for three sample sizes L = 10, 20, 40. 
The distribution is averaged over the disorder, where each 
realization enters the result with the same weight, in- 
dependent of the number of ground-state configurations 
which were available. The distributions extent over large 
intervals down to q — which indicates the existence of 
a complex ground-state landscape. For small overlaps no 
large change is visible with increasing system size, but 
the peak of the distributions located at larger q-values 
shifts to smaller values. 
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FIG. 5. Distribution of overlaps -P(Iq|) for ground states of 
2d ± J Ising spin glass for L = 10, 20, 40. Only for large values 
of 5 a difference is visible, so even for large systems there is 
a finite probability of overlap q = 0. The lines are guides for 
the eyes only. 



1.1 



D" 




0.7 - 
0.6 - 

0.5 I ■ ' ■ ' ■ ' ■ ' 

10 20 30 40 

L 

FIG. 6. Position of the peak of the distribution of over- 
laps -P(|g|) as function of system size for 5 < L < 40. A fit 
function of the form gmax(i) ~ (Jea + qiL~''^ is given. The 
Edwards- Anderson order parameter is obtained by L ^ oo as 
QEA = 0.50(9). 

The position gmax of this peak can be used to calculate 
the Edwards- Anderson order parameter qea which is the 
maximum value of q where P{q) is nonzero in the infinite 
system. Fig. ^ shows the value of (/max as function of 
L. Using a FSS fit with (7max(i) = Iea + <IiL~'^^ a value 
of qEA = 0.50(9) is obtained. The resulting function is 
shown in the figure using a line. 



FIG. 7. Variances a^dgl) of the distributions of overlaps 
Pdgl) as function of L for 5 < L < 40. For L — > oo the value 
converges to a small value. A fit of the form Soo + siL""^ gives 
Soo = 0.004(8). The straight line shows a fit with Soo = 0, 
which should be true for a simple ground-state landscape. 

So far we have seen that finite systems exhibit a com- 
plex ground-state landscape. But to decide whether this 
is true even for the infinite system one must investigate 
the shape of f'd?!) as function of system size L. In fig. 
the variance (see def. (^) of the distribution is shown 
as function of L. By fitting the variance to a function 
of form cr^(L) = Soo + sii"*^ ^ value of Soo = 0.004(8) 
was obtained. This result is very close to zero. In the 
figure a fit with Sqo = is shown which looks very reason- 
able. So it is possible that for L ^ oo the width of the 
distribution shrinks to zero, which would mean that the 
ground-state landscape is simple, as described by the DS 
framework. This is different from the case of the three- 
dimensional model, where the same procedure resulted 
in sf^ = 0.0608(6) [|o). 

However, the impression mediating from fig. |^ is dif- 
ferent: a long tail down to q = persists for all sys- 
tem sizes, so a width of zero for the infinite systems 
seems unlikely. Consider a infinite system, where the 
overlaps are distributed according a distribution with 
a constant probability-density of 0.5 for q < 0.5 and 
a delta-function at q — 0.5 with weight 0.75: P{q) = 
0.259(0.5 - p) + 0.75S{p - 0.5) {p > 0), which seems 
plausible from fig. ^ and ^. Then one obtains a variance 
cr^dgl) — 0.017 which is very close to Soo regarding the 
given error-bar. 

Since the result is not definite so far, next the contribu- 
tion of small overlap- values to P(|g|) is studied. The frac- 
tion Xqg (see def. (Q)) of the distribution below a given 
value qo is displayed in fig. || using qo = 0.2 < 0.5 = qsA 
as function of system size. As Xq^ seems to be indepen- 
dent of the system size it is reasonable to conclude that 
the infinite system has a broad distribution. 
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So far the question whether the ground-state landscape 
is complex has been addressed: a complex landscape for 
the 2d ± J spin glass seems likely, but the results are less 
definite than earlier results for the 3d model. 




In the second part of this section it is studied whether 
the ground states are ultrametrically ordered. For that 
purpose for each realization all possible triplets of ground 
states were chosen and the corresponding three over- 
laps evaluated. The the quantity Sq which is the dif- 
ference between the two smaller overlap values was cal- 
culated (see equation (|ll|)) for all possible triplets with 
a constraint for the largest of the three overlaps: G 
[fgcA - 0.05, ^EA + 0.05] = [0.25,0.35] (which implies 
6q < 0.35 as well). For an ultrametric system the distri- 
bution of Sq should be a delta function. To improve the 
statistics we used the absolute value of all overlaps. The 
distribution P{5q) is shown in fig. ^ for L = 10,20,40. 
Each realization enters the distribution with the same 
weight. One can see that the distribution is independent 
of the system size. The average value of 5q as function 
of system size shown in the inset. Hence, in the infinite 
system > is possible. It seems that the 2d ± J spin 
glass is not ultrametric, which is in strong contrast to 
the result for the three-dimensional model where P{5q) 
converges to a delta- function [EJ . 



FIG. 8. Fraction Xq^ of the distributions of overlaps 
below go = 0.2 as function of L for 5 < L < 40. The value of 
Xo,2 does not decrease with growing L indicating a complex 
ground-state landscape. 
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FIG. 9. Distribution P{Sq) for different system sizes 
L = 10, 20, 40 where Sq = q2 — qi and qi < q2 < f?3 are triplets 
of absolute values of overlaps from independent triplets of 
ground states. Only triplets with qa G [0.25, 0.35] are used. 
For a infinite ultrametric system Sq = holds. The distribu- 
tion turns out to be independent of the system size indicating 
the absence ultrametricity for the ground states. The inset 
shows the average value of Sq as function of system size L. 




-1.0 -0.5 0.0 0.5 1.0 

8q 

FIG. 10. Distribution P2-fix{<l) for different system sizes 
L — 10,20,40 where q G {<?i, 92,^3} and gi < 52 < Is are 
triplets of overlaps from independent triplets of ground states. 
Only g-values of triplets are used where the two other over- 
laps are within the interval [0.25,0.35]. Then for a infinite 
ultrametric system q > 0.25 should hold, while for a metric 
system just q > —0.5 must hold. The small inset shows the 
value integrated from —1 to 0.25. With increasing system size 
the fraction of the distribution below 0.25 remains constant 
although the distribution itself changes. The lines are guides 
for the eyes only. 

Additional information can be obtained by fixing two 
of the three overlaps of a triplet. We took all triplets 
where two arbitrary overlaps fell into the interval [|(7_ea~ 
0.05, \qEA + 0.05]. The resulting distributions J2- fix jq) 
of the third remaining overlap is shown in fig. nffl for 
L — 10,20,40. The triangular inequality gives q > —0.5 
whereas for an ultrametric system q > 0.25 must hold. 
We concentrate on the part of the distribution with q < 
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0.25. Although the shape of the distribution changes 
a httle bit, the fraction of the overlaps forbidden in an 
ultrametric system keeps fairly constant: the inset shows 
the fraction /2_/ia;(0.25) of the distribution below q = 
0.25. Again reasonable evidence against an ultrametric 
organization of the ground-states is found, but it is a 
little bit weaker, since very small values near q = —0.5 
disappear for large systems. 




By studying triplets of ground states with one or two 
of the three overlap values fixed convincing evidences are 
found that the ground states of 2d ± J spin glasses are 
organized not in an ultrametric manner. This is very 
different from the three-dimensional model. To our opin- 
ion this may be connected with the fact that 3d realistic 
spin glasses have an ordered phase at finite temperature 
whereas the 2d model shows a spin-glass phase only at 
T = 0. 

Since only moderate system sizes up to = 40^ were 
investigated, the behavior might change at larger sizes. 
But the finite-size behavior of the data presented here is 
very smooth, thus, we believe that our results persist for 
L oo. 

The results presented here were obtained for a bimodal 
distribution of the interactions. For spin glasses with 
Gaussian distribution the ground state is unique, imply- 
ing that Pidg)!) are both the same in the MF and DS 
picture. The predictions are different at finite tempera- 
ture. We expect the same behavior as for the ± J model if 
one allows deviations of order one from the true ground- 
state energy, but for that system type ground states are 
much harder to calculate using the genetic CEA algo- 
rithm. 
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CONCLUSION 

Many different and independent ground states for 2d 
± J spin glasses were calculated up to sizes L — AO using 
the genetic Cluster-Exact Approximation. From former 
calculations and comparison with exact results is it clear 
that true ground states were obtained. 

By evaluating the distribution of overlaps evidence for 
a complex organization of the ground states is found, but 
the evidence is weak, since it is not clear whether the 
width of the distribution of overlaps scales to zero with 
increasing system size. This is similar to the 3d case, but 
there the evidence for the MF picture is stronger. 
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